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Abstract 

We show the relationship between the mathematical framework used in recent papers by H.C. 
Rosu, S.C. Mancas and P. Chen (2014) [IH3] and the second-order confluent supersymmetric 
quantum mechanics. In addition, we point out several immediate generalizations of the 
approach taken in the latter references. Furthermore, it is shown how to apply the generalized 
scheme to the Dirac and to the Fokker-Planck equation. 
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1 Introduction 

In their recent work M, Rosu et al. present and discuss several interesting examples of 
quantum potentials and their associated zero-energy solutions that are constructed through a 
method inspired by the quantum-mechanical supersymmetry (SUSY) formalism [5]. In general, 
the latter formalism provides a method for the construction of solvable spectral problems that 
are governed by the Schrodinger equation. There is a vast amount of literature on the general 
topic, for a general overview the reader may refer to the reviews [5HZ]- The general SUSY 
formalism splits into two different algorithms, which are usually referred to as the standard 
algorithm and the confluent algorithm. The standard algorithm is a direct generalization of 
the Darboux-Crum transformation w whereas the confluent case can be seen as ingenious 
iterations of the mentioned transformation. Both algorithms serve the same purpose, but they 
are technically inequivalent and as such lead to different results. Despite being less known than 
its standard counterpart, the confluent SUSY algorithm has been subject to studies in a variety 
of articles. The first applications of this algorithm are presented in [IqHIS], while for recent 
progress the reader may refer to [16H22j and references therein. 

The main purpose of this work is to show the close relationship between the mathematical 
method used in m, which for the sake of simplicity we will call “Zero-mode SUSY scheme”, 
and the second-order confluent SUSY algorithm. 

In section [2] we will present the second-order confluent SUSY algorithm and its relationship 
with the zero-mode SUSY scheme, also we will comment on some important features that the 
zero-mode SUSY scheme could take from the confluent SUSY algorithm. The focus of section 
[3] is to apply and extend the results of the previous section to the Schrodinger, Dirac and 
Fokker-Planck equations. Finally, in section [4] our concluding remarks will be presented. 


2 Second-order confluent SUSY algorithm 

In the following we will develop the confluent, second-order SUSY algorithm in order to show 
its close relationship with the zero-mode SUSY scheme, also we will point some important 
properties of the confluent algorithm that can generalize the zero-mode scheme. To this end, let 
us start out with the Schrodinger equation 

^" + (e-Ui)T = 0, (1) 

where the prime stands for the derivative with respect to the position, e for the stationary 
energy and Vi denotes a time independent potential. We will now apply a first-order SUSY- 
transformation to the solution T of ([T|), given by 

— vf 

T =-T -h T', (2) 

u 

where the transformation function u is a solution of 


u" + {\-Vi)u = 0. 


(3) 


This is the same equation as m, but for a stationary energy A (called factorization energy) that 
is allowed to be different from e. Observe further that and u must be linearly independent in 
order to avoid vanishing of ([2]). The function T then solves the Schrodinger equation 


T" + (e - Us) = 0, 


where Us 


Ui -2 


dx'^ 


log (u). 


(4) 


In the next step we apply a first-order SUSY transformation to the new equation Q, which 
requires another transformation function u to generate the function 


^ + (5) 

u 

Let us choose the factorization energy of u to be the same A as for u in ([3|), such that u must 
solve 


u” + (A - Us) u = 0. 


( 6 ) 


Our iteration of the two transformations (I2|), (0) is referred to as confluent SUSY algorithm, 
because the factorization energy is maintained constant [U]. While in the present case we are 
performing a transformation of second-order (two iterations), the confluent SUSY algorithm 
can be applied for arbitrarily high order m- Now, the function T must be a solution of the 
following Schrodinger equation: 


+ (e - Us) T = 0, 


where V 3 


1/2 - 2 ^ log (ii) , 


(7) 


In order to construct the potential V 3 , we recall that a particular solution of equation Q is given 
by Up = Iju. Application of the reduction-of-order formula then yields the general solution of 
([ 6 ]) in the form 


u 


X 

Cl Up + C2 Up ^ dt 

J U^ 


X 



( 8 ) 


where ci and cs are arbitrary constants. Note that in order to obtain ([ 8 ]) we use the fact that the 
differential equation Q does not contain a term with hrst order derivative, in these cases the 


Wronskian of the two linearly independent solutions is constant, and without loss of generality 
we set this constant equal to one since it can be absorbed in the coefficient C 2 - We can now use 
the representation of V 2 in (IH) and the explicit form Q to determine the transformed potential 
V 3 as follows. 


^3 




(9) 


note that every ratio ci/c 2 defines a different potential. Let us now recall that the corresponding 
expression for the transformed potential was obtained in m- Adopting the notation used in 
the latter reference, this potential reads 




dP 



( 10 ) 


where 7 is an arbitrary constant and solves ([T|) for e = 0 (zero-energy solution). If we now 
compare m and m, we observe that Vi-y for 7 = ci/c 2 is a special case of V 3 . While ^0 
can be interpreted as a transformation function of our first SUSY transformation ([2]) for the 
factorization energy A = 0, its counterpart u in Q does not carry any constraints regarding its 
factorization energy. In particular, if our initial Schrodinger equation admits a discrete spectrum, 
\ho represents the ground state solution, while u can stand for the ground state, an excited state, 
or any other solution of (l3|). Let us now construct a solution to the Schrodinger equation ([7]). 
To this end, we must distinguish the cases e = A and e / A. In the first case e = A we know that 
Tp = 1/tt is a particular solution of ([7]). Upon application of the reduction-of-order formula we 
obtain its general solution in the form 




Cl + C2 


^2 

P 


dt 


+ 

u u 


dt 


+ ^ vP dt 


Cl+C2 / 4 -+ 

J \ C2 



( 11 ) 


where Ci and C 2 are arbitrary constants. Since every ratio ci/c 2 defines a potential ([9]) and as 
a consequence a Schrodinger equation ([7|), this parameter is also found in its solutions [13]. We 
will now compare this expression to the corresponding result that was given in mm- It reads 

-• ( 12 ) 

It is straightforward to see that (fT 2 ]) is a special case of ([TT]) if we employ the settings ci/c 2 = 
7, Cl = 1, C 2 = 0 and further identify u = Tq. In the context of spectral problems, m 
is known as the “missing state” solution, because its particular form can be associated with 
a spectral value that has possibly been removed or added by the SUSY transformation. In 






general, all solutions resnlting from a confluent SUSY transformation, can be expressed in terms 
of Wronskians, in particnlar for higher-order cases [ 21 ]. As an example, let us mention the 
Wronskian representation for the solution (|12p : 


^07 




(13) 


where uq and solve the same equation with = 1. The function ui is a solution of the 

inhomogeneous equation 


u" -Vi ui = -Tq. 


Evaluation of (fT^ leads to the identity [7| 


^07 


^0 


^0 

7 + / dt 


which coincides with (| 12 p for a constant of integration 7 . 

In the second case, when e 7 ^ A, the solutions T of equation ([T|) are transformed as follows: 




where ui satisfies the inhomogeneous equation 


(14) 


u'l + (A — Vi) ui = —u. 

The reader might refer to |2n [l211[^ for farther information about the general concept of Wron¬ 
skian representation. Let us now summarize the benefits of embedding the zero-mode SUSY 
scheme into the context of the confluent SUSY algorithm. The principal advantage is that the 
SUSY confluent algorithm provides the possibility of constructing quantum potentials using a 
transformation function u shown in Q that does not need to be the zero-energy mode of the 
initial Schrodinger equation ([T|), but can be any solution of the latter equation. In spectral 
problems, this feature can be used for creating and annihilating bonnd states that lie above the 
ground state level, and also to generate isospectral systems [TOKHlEllESI. Furthermore, it is 
straightforward to extend the scheme used in the latter two references to the case of higher- 
order SUSY transformations, which can be either confluent or a mixture between confluent and 
conventional transformations. Such scenarios have been studied in |21j . 


3 Applications 

We will now present applications of the generalized SUSY scheme that was introduced in section 
2. The focus of our applications will be twofold. We show how our scheme can be used for models 
not governed by a Schrodinger equation, and we choose the transformation function u to solve 
(|3|) for A / 0. In the following applications we will adopt the the notation 7 = ci/c 2 in (|9|) and 
(| 11|1 as a parameter that characterize the transformation. 

3.1 Dirac equation for an inverted harmonic oscillator potential 

In this example we consider the following stationary, two-component Dirac equation in one 
spatial dimension 


i £72 4>' (Vq - E) ^ 


0 


(15) 







where <72 stands for the second Pauli matrix, E denotes a real-valued number, and the spinor 
$ = (<l>i,<I> 2 )^ represents the solution. Furthermore, we assume the potential Vq to be of 
pseudoscalar form, that is. 


Vo = mao + qoai, 


(16) 


for a constant, positive mass m, Pauli matrices cJi, < 73 , and a function qo that is given by 


qo = -x^. 


(17) 


This function can be interpreted as an inverted oscillator potential. In order to construct a 
solution to the Dirac equation (jlSp . we write it in components 

^1 +{E+ m)^2 = 0 , (18) 

$2 - $2 - (-E - m) $1 = 0, (19) 


where we took the form (1161) of the potential and its parametrizing function (1171) into account. 
The system (jlSp . ()19p does not admit closed-form solutions except in the case \E\ = m, which 
we will focus on from now on. We find the general solution of the system to be 


$1 


X 


ki exp 



2 m exp 



J ^2 exp 



dt 


$2 


k2 



( 20 ) 


( 21 ) 


introducing two arbitrary constants ki and k 2 ■ On substituting (ITTP , this solution simplihes to 


$1 = ki exp (—-x^]+-k 2 m exp [ x^ ] x E 2 ( x 

^2 = k 2 exp ( x^ 


( 22 ) 


where E stands for the exponential integral function [26]. We will now apply the SUSY trans¬ 
formation scheme that was introduced in section 2 | 22 ] . Since the latter scheme starts from a 
second-order equation of Schrodinger type, we must first convert the system (fTHI) . (fT9P into such 
an equation. Solving (fTsP for <^2 and substitution of the resulting expression into (fT^ leads to 

$1 + (-x^ + 2 x) d>i = 0. (23) 

Let us note that without incorporation of the explicit form (|17l) . the latter equation would have 
the form 


*^1 + ('?o + do) *^1 — 0- (24) 

This will become important below, when we need to reverse the steps of passing from the Dirac 
equation to its Schrodinger counterpart. The second-order equation (1231) can be compared with 
yielding the identifications 

T = $1 e = 0 Vb = x^ — 2 X, 

recall that $1 is given in (|22p . In order to perform the SUSY transformation (jlll) . we must 
provide a transformation function u that is a solution of equation (1231) . As such, u must have 


the form (j 22 h or it can be any special case of it for a particular choice of parameters fci and /c 2 - 
Due to the length of the expressions involved we omit to state the full form of the transformed 
function (jll[) . Instead, we employ the following parameter setting 


h = 1, Cl = -1/10, 7 = -1/10, k2 = C 2 = 0. (25) 


For these parameters, the function 'h takes the form 

3 exp (—i 


'F = 


3 + 10 X E 2 (I ’ 


(26) 


note that E represents the exponential integral function [26]. The corresponding potential V 3 
that is given in Q reads for the present case 


. 120 exp (—i x^) 1800 exp (—1 x^) 

E3 = x^ - 2 X - J +- ^ V 3 ; 


3 + 10 X E 2 (I x^) 


3 + 10 X E 2 (I x^) 


12 ■ 


(27) 


The solution (|26p and its potential (|27p are displayed in figured] Note that both solution and 
potential are bounded at the origin and nonsingular on the positive real semi-axis. As in the 
case of the function T, the general expression for V 3 is too large to be displayed here. 



Eigure 1: The solution (|26[) in the left plot 
plot. 


E3 



and the corresponding potential (I27p in the right 


Now that we have completed the confluent SUSY scheme, we must relate (f26p and (f27p with 
a transformed Dirac equation that has the form (|15ll . To this end, we first recall that T is a 
solution to ([7|) for the potential V 3 . If we can rewrite the latter potential as 

^3 = qj + q'i, 


then qi will be the parametrizing function of the transformed pseudoscalar potential in a Dirac 
equation of the form (|15p . It is well-known that qi can be constructed by means of the setting 


qi = 


T' 
T ' 


Note that the function T was determined, but we only displayed the particular case (|26p here. 
Eollowing this idea, we show the explicit form of qi for the setting k 2 = 0 


qi 


—x^ -|- 


30 exp (—I x^) 

3 -|- 10 X E 2 (I x^) 


(28) 
















<?1 



Figure 2: The parametrizing function qo (dashed curve), together with its transformed counter¬ 
part Qi for the settings ki = 1, 7 = — 0 . 1 , C 2 = 0 (solid curve). 


Figure [2] shows a particular case of the function qi, together with its initial partner function q^. 

It remains to find the solution of the Dirac equation for the pseudoscalar potential parametrized 
by the function qi, that is, 

i a 2 + {Vi — E) ^ = 0, Vi = m as + qi ai, (29) 

where we assume that \E\ = m. Since the latter Dirac equation has the same form as (1151) . its 
solutions must be of the form ([ 20 l) , (l 2 T]l , where q^ is to be replaced by qi , as shown in (f28]l . Since 
the integrals in the resulting expression cannot be resolved, we omit to show the combination of 
(f28|) and (12(1]) . (|2T]) . Hence, the application of the SUSY scheme to the pair of Dirac equations 
(HSI), (129]) is complete. 

3.2 Schrodinger equation for a constant potential 

We will now pick up an example that was presented in [2] and show that we can generalize it 
considerably when using the full scheme introduced in section 2. Let us consider the Schrodinger 
equation ([T]) for the setting Vi = c^, where c is an arbitrary real number: 

T" + (e - c^) T = 0. (30) 

This equation has the following general solution 

T = fci exp — e + k 2 exp (^—V — e x'j , (31) 

where ki and ^2 stand for arbitrary constants. In the work [2|, the following parameter setting 
was considered 


u = exp (—c x), e = 0, Ci = 1, C 2 = 0. (32) 

Here, the function u is a special case of (1311) for zero energy. As a result of (|32l) . the following 
transformed potential ([9|) and corresponding solution (fTTI) were obtained 


V"3 




2 8 8 

H-1- 

2 c 7 exp (2 c x) - 1 [2 c 7 exp (2 c x) - 1 ]^ ’ 

2 c exp (c x) 

2 c 7 exp (2 c x) — 1 


(33) 


(34) 












This potential and the zero-energy solution coincide exactly with the finding in [ 2 ], except for 
some minor algebraic simplifications. We are now ready to construct a generalized version of 
both pd]) and p^ . While the approach presented in [ 2 ] requires the function u in (l 3 ^ to be 
a zero-energy solution of the initial equation ( 1301 ) . our scheme allows to drop this requirement. 
We will now use the complete solution ( 1 ^ to evaluate ([ 9 ]) and (fTTTi : 


^3 = C^- 




4 (c^ — e) ki exp (^2V— e exp 




27\/— e -|- kf exp ^2\/c^ — e + 4/ci/c2\/c^ — e x — k^ exp 2\/i 
4 (c^ — e) kf exp (^2y/(^ — e + k2 exp (^—2^/c^ — e x^ + 2kik2 


— ex 


27V— e -|- kf exp ^2-\/— e x^ — k^ exp x^ + 4kik2Vc^ — e x 


( 35 ) 


and 


4- = 


2\/c^ — e ki exp ( \/e? — e x] +k2 exp ( —^/e? — e 


— ex 


27 \/— e + k\ exp ^2 \/— e x^ — k"^ exp 2 \/— e x^ + 4/ciA;2\/? 

Cl -|- C2 (^ 2y/ 2 ^2 \/— e x^ — ^2 exp (^—2y/— e — 8kik2{c^ — e) x^ 


-47^ (hk + X + + 

^ ^ kik2 + ^2 exp ( — V — e xj 


, 


( 36 ) 


note that we leave the energy e and the constants ki, k2 undetermined. The last two expressions 
are very general and some particular values for the parameters involved can considerably reduce 
the expressions or be interesting in the physical point of view. 


Potentials with one bound state In equation (1331) can be seen a family of potentials char¬ 
acterized by the parameter 7 and (f34]) presents its only bound state at zero energy. Now let us 
consider the following setting 

u = exp — e , e < (? ,Ci = 1 , (72 = 0 . ( 37 ) 


The potentials obtained after the confluent SUSY transformation and its bound state are 

8(c2 - e) 8(c 2 - e) 


U3 = c^- 


4^ = 


27 ^^ eexp(—2\/c^ — e x) -|- 1 

2 y/— eexp(\/c^ — e x) 

2^y/(? — e -|- exp(2\/ e? — e x) 


-b 


27 ^^ eexp(—2\/c^ — e x) -|- 1 


,( 38 ) 


( 39 ) 


In order to avoid singularities in the new potential 7 has to be a positive constant. The function 
(fMl) is the only bound state of the system. Figure [ 3 ] shows the graph of the potential Vi = e? 
(dotted curve) and two potentials V3 with different parameters. In the first case we use a 
zero-mode solution as transformation function (dashed curve) and then a solution of ([ 3 ]) where 
A = 0.6 (solid curve). In both cases the systems have a bound state associated to e = 0 and 
e = 0.6, respectively. 




























































Figure 3: The constant potential Vi = 1 (dotted curve) and the potentials V 3 with the parameters 
e = 0 , 7 = 200 (dashed curve) and e = 0 . 6 , 7 = 200 (solid curve). 


Hyperbolic case Simplihcation of (j35p and (I36p considering ki = k 2 = Ci = 1, C 2 = 0 and 
also e < (? gives the following results 


2 16(c^ — e) cosh(\/— e x) sinh(\/— e x) 

'J(? — e (7 + 2 x) + sinh( 2 V (? — e x) 


T = - — 

\/(? — e(7 + 2x) + sinh(2Vc^ — e x) ’ 


16(c^ — e) cosm(Vc^ — e x) 

\l(? — e (7 + 2x) + sinh(2\/— e x) 
2\/— e cosh(\/— e x) 


(40) 

(41) 


where for the sake of simplicity we introduced hyperbolic functions. The family of potentials is 
singular for all the values of 7 and as a consequence no bound states can be found. 


Trigonometric case When the parameter e is greater than the potential Vi = (? the exponen¬ 
tial functions in (j35l) and (I36p can be expressed as trigonometric functions. With the following 
setting 

u = sin — e x^ , e > c^, Ci = 1, C 2 = 0, (42) 

the transformed potential and the corresponding solution of the Schrodinger equation are 


2 8(c^ — e) sin(2-\/c^ — e x) 

2\/(? — e (27 -|- x) — sin( 2 \/c 2 — e x) 

4 \/(? — e sin(\/c^ — e x) 

2\/(? — e (27 -|- x) — sin( 2 Vc^ — e x) 


32(c^ — e) sin^(-\/c^ — e ic) 
2 \/c^ — e (27 -|- x) — sin( 2 -\/— e x) 


(43) 

(44) 


If we consider the whole real axis as domain then there are no regular potentials V 3 obtained 
with the confluent SUSY transformation. 


3.3 Fokker-Planck equation with constant diffusion coefficient 

The SUSY scheme introduced in section 2 can also be applied to certain cases of the Fokker- 
Planck equation, as will be shown in this subsection |27ll28j . Let us consider the following 







































equation 


^ e (0,oo) X (0,oo). (45) 

This Fokker-Planck equation (|45p has a constant diffusion coefficient and a drift function U that 
we assume to depend on the spatial variable only. We will comment on boundary and initial 
conditions once we have constructed the transformed equation and its solutions. In order to 
apply our SUSY algorithm, we must convert (f45|) into a Schrodinger equation of the form (fT|). 
To this end, we rewrite the solution / as follows: 

/ = exp {—U — k t) 'k, (46) 


where /c is a real constant and T is a function independent of the time variable. Substitution of 
(1461) into equation (|45p and simplification gives the following result 

T" +[2k- {U'f + U''] T = 0, (47) 

This matches our Schrodinger equation ([1]) if the following settings are employed 

e = 2k Ui = {U'f - U". 


In the next step we specify the drift function U in (1451) as follows: 


" - T- 


Let us now rewrite equation ()T7p using the explicit form (|4SP of U : 

T" + [2 fe - + l] T = 0. 

The general solution of this equation is given by 


(48) 


(49) 


(50) 


where ki , k 2 are free constants and H stands for the Hermite function [26] . Recall that for 
an index that is a nonnegative integer, the Hermite function becomes a polynomial. In the 
next step, we must choose a solution u of (|49p that enters in both the potential dOj) and the 
corresponding function (fTH . We pick the following special case of (1501) 


T = fci exp I- Hk{x) + k 2 exp I — H_k-i{i x 


X 


u = exp 


Hk{x), 


(51) 


where k is a nonnegative integer. Substitution of (f5T|) into ([9|) yields 

2 exp (-2 x^) H^{x) 


U 3 = x"-l + 


1 2 


+ 


7 + /exp(-t 2 ) Hl{t) dt 

4 exp (—x^) X H‘1{x) — 8k exp (—x^) Hk-i{x) Hk{x) 


(52) 


7 + /exp(-t 2 ) Hl{t) dt 

The solution of our Schrodinger equation ([7|) can be found by inserting ([5X1) into ([TT]) . We get 


4' = 


exp 


Hk{x) 


7 + f evip {-U) Hl{t) dt 


C1 + C2 


exp (-t^) 


-\ 2 


7 + / exp (-s^) iL|(s) ds 


dt 













( 53 ) 


In the last step we need to convert (|53p into the solution of a Fokker-Planck equation, the drift 
function of which is determined by (1521) . To this end, let us introduce the following quantity 

F = -log('k), (54) 

where 'k is defined in (1531) . It is then straightforward to verify that the Schrodinger equation 
([7|) for the potential (l52l) can be rewritten as follows 

4>" +[-{V'f + V'']^ = 0. (55) 

Note that for the sake of clarity we did not include the full form of V, which can be obtained 
from combining (1531) and (|54p . Now, comparison of (1551) with (|47l) with (|54ll shows that the 
function 


g = exp (—I/)'!' = 


(56) 


is a time-independent solution of the Fokker-Planck equation 


dg 1 d'^g 
dt 2 dx^ 



0 . 


(57) 


We observe that (j54j) is the transformed drift function, while the corresponding solution (I56p is 
obtained by squaring 'k in (j53p . Since the integrals in the latter expression cannot be resolved 
in their general form, we will now restrict ourselves to a particular case. Let us employ the 
settings 


7 = -0.9, Cl = - 0.25 (72 = 0 


This gives the transformed drift function as 


V 


-log 


5 exp 

18 — 10 yTr erf(x) ’ 


k = 0. 


(58) 


(59) 


where the symbol erf stands for the error function [2^. The left plot in figure 0] shows the graphs 
of the drift functions V and U, as given in (f59P and (f48p . respectively. 

Inspection of expression (I59p shows that the transformed drift function V can be defined on 
the whole real line. The same holds for the associated solution of the transformed Fokker-Planck 
equation ([57P , see the right plot in figure 01 In addition, we observe that this function vanishes 
at the infinities.A study of expression (j54j) indicates that for 7 S (— 00 ,0) when A: = 0 and Ci < 0 
will lead us to a regular drift function V. 

Let us now choose a set of parameters that is different from ()58p 

7 = 0 . 1 , Cl = 1, C 2 = 0, k = 1 . (60) 


For these settings, the transformed drift function reads 


V 


-log 


20 X exp 


—20 X + exp (x^) [1 -|- 10 erf(a:)] 


(61) 


Similar to the previous case, the left and right plots in figure [5] show the drift functions (08]), 
dsm and the associated solutions of (05|), ([57p . respectively. We observe that the transformed 
solutions can be defined on the positive half-axis, vanishing at zero and at positive inhnity. 
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Figure 4: Left plot: the transformed drift function (|59p (solid curve), together with its initial 
counterpart (f48|) (dashed curve) for the settings ([58]) • Right plot: solutions of the Fokker-Planck 
equations ([571) (solid curve) for t = 0 and ([1^ (dashed curve). 




Figure 5: Left plot: the transformed drift function ()6ip (solid curve), together with its initial 
counterpart (I48p (dashed curve) for the settings (|6UI) . Right plot: solutions of the Fokker-Planck 
equations ([571) (solid curve) for t = 0 and ([1^ (dashed curve). 


Furthermore, a study of equation (l54p under the set of parameters Ci = 1, 6*2 = 0, k = 1 shows 
that when 7 S (0, 00 ) the only singularity of the drift function V is in the origin. 

Before we conclude this example, let us show some further solutions in hgure [6] that fulfill 
boundary conditions on the whole real axis (left plot) and the positive half-axis (right plot). 

We observe that these solutions correspond to excited states of the Schrodinger equation ([7]), 
which are not accessible through the method introduced in m, but only after its generalization 
presented in section 2. 

4 Concluding remarks 

We have shown the close relationship existing between the zero-mode supersymmetry scheme 
and the confluent SUSY algorithm of second order. With this relationship we were able to take 
advantage of some features of the confluent algorithm in order to generalize the results of H. 
Rosu et. al. (2014) [T]-[3]. One of these features is the fact that there is no need to use only the 
zero-mode solution in order to perform the transformation. Moreover, we applied the second- 
order SUSY confluent algorithm to three different equations. First, we took a one dimensional 
Dirac equation for an inverted harmonic oscillator potential and performed a confluent SUSY 
transformation in order to generate a new Dirac equation. As a second example a Schrodinger 


















Figure 6: Solutions of the Fokker-Planck equations (1571) (solid curve) for t = 0 and (|15]) (dashed 
curve) for the settings Ci = 20, C 2 = 0, k = 2, 7 = 20 (left plot) and Ci = 20, C 2 = 0, k = 
3, 7 = 100 (right plot). 


equation with constant potential was taken as an initial potential to generate new ones using no 
zero-mode transformation functions. Finally, we showed how to apply the second order confluent 
algorithm to the Fokker-Planck equation using not only the ground state but also excited states 
as transformation functions to generate new Fokker-Planck equations and its solutions. 
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